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Proper motion of 7-rays from microhalo sources 
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I discuss the prospects of detecting the smallest dark matter bound structures present in the 
Milky Way by searching for the proper motion of 7-ray sources in the upcoming GLAST all sky 
map. I show that for dark matter particle candidates that couple to photons the detection of at 
least one 7-ray microhalo source with proper motion places a constraint on the couplings and mass 
of the dark matter particle. For SUSY dark matter, proper motion detection implies that the mass 
of the particle is less than 500 GeV and the kinetic decoupling temperature is in the range of [4-100] 
MeV. 

PACS numbers: 95.35.-l-d,98.80.-k, 95.55.Ka, 98.70.Rz, 97.10.Wn 



In the ACDM cosmological model, ^ 26% of the 
present-day energy density of the Universe is in the form 
of dark matter which remains an outstanding prob- 
lem in cosmology and particle physics. Favored Cold 
dark matter (CDM) particle candidates (see 0| and ref- 
erences therein) are Weakly Interacting Massive Particles 
(WIMPs). They are kept in thermal equilibrium through 
interactions of the form XX ^ ffi ^-nd in kinetic equilib- 
rium through xf ^ xf- The temperature at which these 
particles decouple from the radiation fluid is called the 
freeze-out temperature, while the temperature at which 
they decouple from scattering interactions is called the 
kinetic decoupling temperature (Td) ^\^\^ 

In the CDM paradigm, the first collapsed and viri- 
alizcd dark matter structures, microhalos, are formed 
at high redshiftsj^ 0, @]- Their minimum mass is set 
by the RMS dark matter particle velocities set at ki- 
netic decoupling; scales smaller than the free-streaming 
scale are smoothed out. For the case of supersymmet- 
ric (SUSY) cold dark matter (CDM), this mass scale is 
Mm « IO-^Mq (rrf/10MeV)-3 Recent numeri- 

cal simulations of the formation of microhalos in the con- 
text of SUSY CDM provided hints on their properties, 
and on their evolution and survival in the hierarchical 
structure formation scenario [l(l llll | . Their potential de- 
tection using 7-rays has been studied in Refs. 0, 0| . 
Based on simple analytical studies, it has been argued 
that the abundance of microhalos in the Milky Way can 
be severely reduced due to interactions in early rapid 
merging processes and in stellar and galactic disc inter- 
actions [13, llJ, UM, UM ■ These arguments are still under 
debate and need to be resolved through detailed numer- 
ical simulation studies. 

Investigating the existence of microhalos in the Milky 
Way could be very important for two reasons: 1) A possi- 
ble detection can provide information about the particle 
physics properties of the dark matter particle, and 2) a 
measured abundance in the Milky Way halo contains in- 



formation on the hierarchical assembly of dark matter 
halos at very early times, a task that is unattainable by 
any other method. If a fraction of microhalos survive 
in the present day Milky Way halo, and the dark matter 
particle couples to photons (as is the case for SUSY CDM 
as well as candidates that originate from theories with 
extra-dimensions), they may appear as 7-ray sources in 
the sky with potentially detectable proper motion ■ A 
search for proper motion of unidentified 7-ray sources can 
provide important information on the physics of the dark 
matter particle and on the survival rate of microhalos in 
the Solar neighborhood. In this paper, I present proper 
motion signals that should be expected in the upcoming 
GLAST 7-ray experiment | http://glast.gsfc.nasa.gov I for 
dark matter particle candidates that couple to photons. 
I find that for SUSY CDM particles with mass less than 
< 500 GeV and kinetic decoupling temperature of 
Td ~ [4 - lOOjMeV there wiU be at least one detectable 
microhalo with proper motion if 20% of microhalos at 
the Solar radius survive 0, 0| . 

Photons from Microhalos: The number of photons with 
energy greater than iJth on Earth from a microhalo at 
distance V is 
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Here, {av) is the annihilation cross section of the dark 
matter particle to photons, is the mass of the dark 
matter particle and = J^^[dN~^/dE]dE is the to- 
tal number of photons emitted above a threshold energy 
Eth and dN-y/dE is the spectrum of the emitted pho- 
tons. The effective area of the detector is ^eff, and Texp 
is the exposure time. The term g[p{r)] contains all the 
information about the distribution of dark matter in a 
microhalo at a distance P, g[p{r)] = T)^^ J p^{r)d^r — 

2tt J^"^ smip{J p^[r(s, ip, V)] ds} d-ip. The integral is per- 
formed along the line of sight, for an angular extent of 
4>rn = tan~^[rs/'D], where Ts is the scale radius. The dis- 
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tribution of mass is p(r) and is assumed to be described 
by an NFW profile [13 (for an NFW profile, -90% of 
the flux comes within the region inside of r^). All mi- 
crohalos are described by the same value of the quantity 
o[p(r)] if all are formed at the same redshift, Zform 
|lOj. Given a concentration parameter Cv ~ 1 W], a 
microhalo of mass is described by an NFW pro- 
file with a characteristic density of ps ~ 756 GeV cm~^ 
and a scale radius « 3 x 10^^ pcM^^'^^g^^'^, where 
Af,„,_6 = M^/IO-'^Mq and Z70 - [(1 + Zform)/71] For 
the case of SUSY CDM, the number of photons detected 
on Earth that originate within a microhalo region defined 
by the PSF of GLAST can be written as 
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Here, Pq.os = 2?/0.05pc, the threshold energy is Eth = 
3 GeV and the Hubble parameter is taken as h = 0.7. In 
this expression, {av)-2e = (CTi;)/10~2^cm'^s~^, = 
Mj^/40 GeV, and the photon spectrum is taken from 
Ref. 01 . This particular choice of cross section, particle 
mass and threshold energy represents an optimistic ap- 
proach to the case of SUSY CDM detected with GLAST 
[l9j, because the annihilation cross section can be up to 
a factor of 10^ smaller, while the mass can be as high 
as numerous TeV. However, for dark matter candidates 
that arise from stable Kaluza-Klein excitations in mod- 
els with universal extra dimensions the annihilation cross 
section to photons can be much higher 0]. 

In order to detect the 7-ray signal from microhalos, the 
source signal must be larger than the threshold number of 
photons for detection, i.e., N.y^s > N^^th, where iV^.th ~ 
6 for a GLAST orbit-averaged effective area of AcS ~ 
2 X lO'^cm^ [2^ and an exposure time of 2 years. Given 
a microhalo of mass Mm, the maximum distance from 
the Sun "Dmax at which a microhalo can be detected is 
obtained by setting Eq. ^ equal to N^^th- I define the 
"visibility volume" of a microhalo of mass Mm to be the 
volume defined through I?max- 

Microhalos formed at Zform ~ 70 will merge to form 
present day dark matter halos. A fraction of the Milky 
Way halo mass perhaps remains in the form of micro- 
halos formed at very high redshifts. Simulations show 
that the subhalo mass function in Milky Way size ha- 
los is dN/dM ~ M~2, with about 10% of the mass of 
the Milky Way in halos of mass greater than IO^Mq. 
Analytical arguments based on approximations for the 
survival pro bability of microhalos in an evolving par- 
ent halo |15|. as well as preliminary numerical work on 
the mass function of subhalos down to microhalo scales 
[ll| . find the form of the subhalo mass function is pre- 
served down to microhalo scales. Thus, the amount of 
mass per logarithmic mass interval is a constant, and 
depends only on the survival probability of microhalos. 
If approximately 20% of microhalos survive jl5l | in the 



Solar neighborhood (with perhaps a small mass depen- 
dence on microhalo mass |l6|), then the fraction of the 
local mass in microhalos is roughly of order 0.2%. I de- 
fine the quantity ^ to be the fraction of the dark matter 
density in the solar neighborhood that is in the form 
of microhalos, i.e., f — 0.002. If there are micro- 
halos in a logarithmic interval in Mm, the number of 
microhalos per radial distance interval from the Sun is 
dNtn/dT> = Att^V^ p{rQ) /Alyn, and the total number of 
microhalos per logarithmic mass interval in Afm within a 
visibility volume is TVm = /J'""" [dA'^m/c^^^Jdl). The Sun 
is located at « 8 kpc from the galactic center, andthe 
local dark matter density is p{rQ) — 0.01 Mq pc^"^ 2l]. 

Microhalo proper motion: I assume microhalos obey 
a Maxwell-Boltzmann velocity distribution function set 
by the virial temperature of the Milky Way halo. On 
the Earth's frame, the distribution function can be writ- 
ten as f{v)dv = [v(^{v)/iT^^'^VeVo]dv. where, (^{v) — 
exp[{vi/vo)^]-exp[{v2/vo)'^], vi = v-Ve and 1^2 = v + Ve- 
Here, vq « 220 km s^^ is the velocity of the Sun in the 
frame of the Galaxy, and Ve — 1.12i;o is the speed of the 
Earth (the very small annual modulation of the Earth's 
velocity is neglected). A microhalo at a distance V from 
the Sun, and with a velocity with magnitude v at an angle 
9 relative to the Earth's frame will exhibit a proper mo- 
tion of angular displacement p, = tan~^[Tv sinO/T)] in a 
time period r. Naively, for GLAST, a detection of proper 
motion requires that the angular displacement must be 
greater than the point spread function, i.e., p, > no = ^ 
arcminutes (which is the 68% PSF for GLAST). However, 
the PSF of GLAST improves to less than 9 arcminutes 
at energies greater than few GeV because the sensitiv- 
ity of the instrument is not background dominated. It is 
reasonable to assume that a detection of proper motion 
is defined by the condition pth > Mo/ \/^'f,th ~ 4 arcmin. 
This is a conservative approach, because large number of 
photon counts (as might be the case for a nearby micro- 
halo) ma y in fact increase the angular resolution of the 
detector |22l |. Furthermore, it is possible that some mi- 
crohalos are close enough that they will have an angular 
extent (see Fig. ^ such that the excess count in adja- 
cent angular bins may be used to even better localize the 
source, increasing the sensitivity of the detector to proper 
motion measurements. The number of microhalos in a 
logarithmic interval in Mm detected with GLAST with 
proper motion greater than pth is obtained by integrat- 
ing the velocity distribution function with the abundance 
of microhalos with a velocity vector at an angle greater 
than OminipthjV,!)) = sin^"^[I?tan/^th/'yT], as 



A^m(Mm,> Mth) = 47r 



(3) 



dV 
X (^J f{v)cos[0niin{pth,v,V)]dv^ dV 

Results: Fig. ^ shows the number of visible microhalos 
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FIG. 1; Top: The dependence of the number of visible mi- 
crohalos per logarithmic mass interval to the microhalo mass 
and the properties of the dark matter particle. Sohd lines 
depict the iso-number contours of microhalos which also ex- 
hibit detectable proper motion. The dashed line shows the 
best case scenario for SUSY CDM particle. The shaded area 
shows the region that is already being excluded by EGRET 
measurements. 



per logarithmic mass interval and contours of the num- 
ber of microhalos with detectable proper motion for dif- 
ferent microhalo masses and physics of the dark matter 
particle. The dashed line depicts the best case scenario 
for SUSY CDM, namely, a neutralino dark matter, with 
= 46 GeV, and (av) = 5 x IQ-^^cm^s-^ The photon 
spectrum is taken from Ref. and is integrated from 
ii^th = 4 GeV. The integration time is 2 years. If, for 
example, the decoupling temperature of the SUSY CDM 
particle is ^ lOMeV, then microhalos have a mass of 
~ 10-" Mq, and if - 0.2% of them survive in the Milky 
Way halo to the present day, ~ 20 will be detected with 
GLAST j an d 8 out of 20 will exhibit a detectable proper 
motion |27l |. 

The shape of the contours of microhalos with de- 
tectable proper motion can be understood with the help 
of Fig. El The radius from the Sun at which a microhalo 

1 /2 

is detected is I? ~ A/m . However, the radius of the vol- 
ume that is needed in order to encompass at least one 

1 /3 

microhalo is I? ~ Mm (see e.g. the intersection of the 
A^(< V) = 1 contour and the detectability line in Fig.[21). 
Therefore, for small mass microhalos < 10-^Mq, the 
number of microhalos with proper motion is limited by 
the non-abundance of microhalos within a visibility vol- 
ume. On the other hand, for high mass microhalos, even 
though the detectable number can be high, the number 
of microhalos exhibiting proper motion is limited by the 
distance at which a microhalo can be placed and still ex- 




FIG. 2: The angular size of microhalos as a function of mi- 
crohalo mass and distance for the best case scenario of SUSY 
CDM. The solid line defines the distance threshold for de- 
tection. Dashed lines depict the iso-number contours of the 
number of microhalos per logarithmic mass interval within a 
volume defined by D. Dot-dashed lines show the maximum 
distance at which a microhalo could exhibit proper motion of 
1, 3, 6 & 9 arcmin. The cross-hatched area corresponds to the 
region where microhalos have an angular extent that should 
be detectable by GLAST. 



hibit a detectable proper motion, because /i ^ D^^ (see 
e.g. the intersection of the fi = 3' and A^(< 2?) = 1 
contom- in Fig.|51). 

For smaller (larger) values of ^, such as the case where 
microhalos are disrupted more(less) efficiently, the num- 
ber of visible microhalos and/or ones that exhibit de- 
tectable proper motion should be multiplied by a factor 
of ^/0.002. An interesting outcome of the detection of 
a large number of microhalos will be the potential mea- 
surement of their parallax. Such measurements can be 
used to estimate the absolute luminosity and thus the 
measured flux can provide important information about 
the physics of the dark matter particle and its distribu- 
tion. This calculation assumes isotropy of the microhalo 
distribution in the solar neighborhood and thus isotropy 
in the distribution of detected microhalos. However, the 
Doppler effect due to the motion of the Earth-Sun system 
in the Milky Way halo will create a dipolar distribution 
with an amplitude of order AN^/N^ ~ 10"^% |2^ . 

In the context of SUSY CDM, the detection of at least 
1 microhalo with proper motion has an important impli- 
cation: the decoupling temperature must be in the range 
~ [4 — 100] GeV, and the mass of the particle must be 
< 500 GeV. The former is inferred from the crossing 
of the SUSY Nj{av)/M'^ line with the > 1 iso-number 
contour of detectable microhalos with proper motion in 
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FIG. 3: The dependence of the number of visible microhalos 
per logarithmic mass interval on the microhalo mass and red- 
shift of formation. Solid lines as in figure Q The dashed line 
depicts the redshift at which microhalos form in numerical 
simulations [To l ln l lT ^ . 



Fig. n] The latter is inferred from the lowest value of 
the quantity N^{av)/M^ along the same > 1 iso-number 
contour in Fig. ^ These bounds can be understood in 
the following way: High decoupling temperatures imply 
small microhalos, which are less luminous; The visibility 
volume shrinks reducing the number of detectable mi- 
crohalos {N^ ~ M^^V^ ^ Mm ^). On the other hand, 
low decoupling temperature results in more massive mi- 
crohalos, in a larger visibility volume. Proper motion is 
inversely proportional to the distance of the microhalo, 
fi 'D~^ ~ Mm^''^, and so more massive microhalos 
will tend to exhibit smaller proper motions. It will be 
of interest to deduce the implications of the kinetic de- 
coupling temperature and particle mass constraints to 
direct searches for dark matter, because the interactions 
(scalar scattering) that set the kinetic decoupling tem- 
perature are the same as the ones that are expected to 
take place in a direct dark matter detection experiment. 

The cross-hatched region in Fig. [21 depicts the region 
where microhalos can potentially be resolved as extended 
objects (greater than the 9 arcmin PSF of GLAST). If 
for example microhalos have mass of order ~ 10^^ Mq, 
then there will be ~ 25 microhalos visible, out of which 
~ 7 will be resolved as extended objects, and ^ 2 of 
them will exhibit a detectable proper motion of greater 
than 3 arcmin. This number is a conservative estimate 
since the excess photon flux from adjacent resolution bins 
in extended objects could be used to better localize the 
position of the source [isj. 

Fig. 13 shows the number of visible microhalos and 



contours of the number of microhalos with detectable 
proper motion for different microhalo masses and for- 
mation redshift. Recent simulations suggest that micro- 
halos are formed at ^form ~ 70 (depicted as the dashed 
line in Fig. O. The assumption that microhalos form 
at Zform ^ 70 should be taken with caution because the 
redshift of formation depends on 1) the rarity of the over- 
density from which the microhalos collapse 

and 2) 

the scale invariance of the primordial power spectrum 
Furthermore, simulations of small box-sizes are 
very sensitive to the redshift at which the simulation is 
initiated [2^. This redshift must be high enough such 
that the initial Zel'dovich displacement is a small enough 
fraction of the mean inter-particle separation. For sim- 
ulations aiming at the formation of microhalos this cri- 
terion implies a starting redshift much higher than what 
has been assumed in refs. [lol Not fuUfiling this 
requirement delays structure formation on all scales in 
simulations. Therefore, if a population of microhalos is 
formed at a higher redshift than what has been assumed 
in this calculation would imply that the number of vis- 
ible microhalos and microhalos with detectable proper 
motion could be higher as shown in Fig. |2] 

In summary, a detection of the proper motion of 7- 
ray sources would be of profound importance. First 
and foremost, it will be a detection of annihilating dark 
matter. In addition, it will provide information on the 
survival rate of microhalos in the Solar neighborhood 
US m El El El El- Furthermore, it wiU provide an 
upper bound to the dark matter particle mass, namely, 
< 500 GeV, an area in the parameter space that will 
be accessible for further studies at the upcoming LHC. 
Moreover, the kinetic decoupling temperature of the dark 
matter particle must be in the range — [4— 100] MeV, 
dramatically restricting the current allowed range of val- 
ues ,26j . Therefore, it is essential that a particular analy- 
sis for the proper motion of 7-ray sources is performed by 
the GLAST team, and that the GLAST lifetime is max- 
imized so as to increase the baseline for proper motion 
detection. 
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